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Abstract 

The problem of synchronization of coupled self-oscillators by ex- 
ternal force is studied. The charts of Lyapunov's exponents in the 
"frequency - amplitude" parameter plane are obtained within the 
framework of the phase approximation. We identified the character- 
istic configurations of the domains of complete synchronization, two- 
and three-frequency quasi-periodic oscillations and different variants 
of partial synchronization of oscillators by an external force. The dif- 
ferences between regimes of mode locking and beats of partial oscilla- 
tors are revealed and discussed. To visualize and analyze the domains 
of resonant two- frequency tori we construct the charts, in which re- 
gions of dynamics with different winding numbers are represented by 
colors. 
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1 Introduction 

In nonlinear science, there is a fundamental problem concerning creation and 
evolution of multi-frequency quasi-periodic oscillations. Such phenomena 
are common in nature and technology, including electronics, laser physics, 
mechanics, as well as neurodynamics, biochemistry, climatology [1]-|5]. 
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In due time, a scenario of onset of hydrodynamic turbulence was advanced 
by Landau and Hopf [B] , [7] based on the idea of sequential appearance of in- 
creasing number of oscillatory components with incommensurate frequencies 
due to successive bifurcation in the course of increase of a control parameter 
(the Reynolds number). Latter, it was criticized and reconsidered by Ruelle 
and Takens [S]. They asserted that after a birth of a few oscillatory modes 
the motion will be associated with a strange attractor characterized by sen- 
sitivity of orbits to initial conditions and continuous spectrum. Although 
the Ruelle - Takens motivation was based on rigorously proven mathemati- 
cal theorem, its application to concrete systems with complex dynamics in 
physics, technology, and other disciplines remains questionable. Particularly, 
many researchers outline typical occurrence of multi-frequency oscillations in 
various examples of nonlinear systems One must acknowledge that 

at the moment we do not have a complete and satisfactory general picture of 
emergence and disappearance of multi-frequency quasiperiodic motions, the 
role of synchronization and chaos in this picture, illustrative visualization of 
respective phenomena etc. 

A convenient class of models for discussion of many aspects of the prob- 
lem is represented by ensembles composed of oscillators (autonomous or non- 
autonomous). It is remarkable that novel and unexpected features of dynam- 
ical behavior of such models were found recently even in the simplest case 
of two coupled self-oscillatory elements, say, when the degree of excess over 
the thresholds of Andronov-Hopf bifurcations are essentially different (see, 
e.g. [18]-[20]). The problem of the dynamical behavior of three interacting 
oscillators is considerably more complex and multifaceted. It can be studied 
in the framework of different approaches: in terms of the original differ- 
ential equations of the system, on the level of slow-amplitude description, 
and with the help of phase equations and model maps for the phases (the 
maps on a torus). For example, in Ref. |2T] the model is discussed com- 
posed of coupled maps of rotations. Within its framework a large number 
of possible types of bifurcations were revealed, as well as the possibility of 
occurrence of "toroidal" chaos was outlined. Note that the problem of the 
dynamics of phases on the torus is complex and deep; in particular, such 
thoughtful and non-trivial approach like the renormalization group method 
was developed in this context ^22], [23]. Theoretical and experimental stud- 
ies of three-frequency quasiperiodic motions in coupled electronic oscillators 
were carried out e.g. in Refs. [I2]-[T3]. One result is that emergence of three- 
frequency oscillations is associated often with the bifurcations of merging and 
disappearance of pairs of stable and unstable two-frequency tori. In Ref. [U] 
the authors investigated a model of two van der Pol oscillators with reactive 
coupling excited by an external periodic force. For weak coupling the dom- 
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inant three-frequency quasiperiodicity was observed, and for large couplings 
chaos becomes possible and typical. Dynamics of a ring of three coupled 
phase oscillators was discussed e.g. in Refs. [2l]-[26]. A possibility of three- 
frequency quasiperiodic motions was noticed in the model of three-coupled 
Lorenz systems [T6], [17], and its destruction accompanied with transitions 
to chaos was considered. 

In recent years, the problem of forced synchronization of two oscillators 
by an external force has been studied in some details. In a series of papers 
[27j-[3T] mechanisms for synchronization of resonance cycles on a torus in 
the autonomous systems were established and discussed. In [30j, |31j phase 
equations describing the excitation of two coupled dissipative oscillators by 
periodic driving force were obtained and analyzed. A picture of domains of 
distinct dynamical behavior on the parameter plane of the frequency and 
amplitude was revealed; it contains e.g. regimes of complete locking of the 
oscillators by the external force, and regimes of two- frequency and three- 
frequency quasiperiodic motions. This result, to some extent, generalizes the 
concept of Arnold's tongues on the situations of multi-frequency motions. 

In the present paper, we study a model of periodically driven system 
composed of two coupled van der Pol oscillators, and corresponding phase 
equations. In comparison with ^27]- [31], one novel respect concerns consider- 
ation of regimes associated not only with the mode locking, but as well with 
beats of the partial oscillators subjected the action of the periodic driving 
force. Bifurcation analysis of the complete synchronization state is carried 
out on the base of relatively simple model describing phase dynamics of the 
coupled oscillators. We perform analysis of the spectrum of Lyapunov's ex- 
ponents in the parameter plane of frequency mismatch and amplitude of the 
external force, and summarize it by means of visualization with the so-called 
"charts of tori" . On such a chart, the resonant two-frequency tori with differ- 
ent winding numbers, which occur on the surface of a three-frequency torus, 
are shown in colors. With this technique, we detect arrangement of the syn- 
chronization domains in the parameter plane ignificantly different for two 
situations: when the isolated two oscillators are mode locked, and when they 
demonstrate the beats. Finally, we compare and discuss typical Lyapunov 
charts in the parameter planes plotted for the approximate phase equations 
and for the original differential equations. 
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2 Equations for the periodically driven phase 
oscillators 



Let us consider a system of equations describing dynamics of two dissipative 
coupled van der Pol oscillators under external harmonic force acting on the 
first oscillator: 

X — (A — x'^)x + (1 )x + — y) = B sin(a;t), 

2 A (1) 
y-{X- y^)y + (1 + — )?/ + ^l{y - x) = 0. 

Here A is the parameter responsible for the Andronov - Hopf bifurcation 
in the uncoupled oscillators; A is the detuning parameter, proportional to 
a difference of the oscillator frequencies; fi is the coupling parameter; B is 
the amplitude of the external force and a; = 1 + is its frequency. Central 
(mean) oscillator frequency is assumed to be normalized to unity. Hence, 
Q is the dimensionless detuning of the frequency of external force from the 
central frequency. 

We use the method of slow amplitudes and present the solution for 
dynamical variables x and y in the form: 

X = ae'^' + a*e-^'^*, y = ce''^* + 0*6"^'^* (2) 

where a{t) and c{t) are complex amplitudes of the oscillators and u is the 
external forcing frequency. Further we use the additional conditions standard 
for the method of slow amplitudes: 

ae*"^* + d*e-''^' = 0, ce*'^* + 0*6"^'^* = 0. (3) 

Substitution of ([2]) and ([3]) into Eqs.(IT]), with multiplication by e"""^* and 
subsequent time averaging leads to a set of truncated equations: 

A B 

2d = Aa — \a\ a — 2i(Q H )a — n(a — c) , 

. , ^A ^ (4) 

26 = Ac — |c| c — 2i{Q — — )c — /i(c — a). 

In the last equation we took into account that a; = 1. Parameter A may be 
eliminated by renormalization. 

Further let us set a = Re^"^^ and c = re*'^^, where i?, r are real amplitudes 
and ^/'i,'02 are real phases of the oscillators determined relatively to the 
external signal. Following [1] , we assume in the last equation that oscillators 
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are moving in the neighborhood of the hmit cycle R = r = 1. Finally we 
obtain the equations: 



i>i = --r - 1^ + -sin(?/;2 - + 6sin(?/'i), 
^2 = J - + - sin(?/'i - -02), 

where ^ = These are phase equations for the periodically driven coupled 
oscillators EU- 



In the case of uncoupled oscillators (/i = 0) the system ([5]) splits into 
two independent equations: phase equation for the periodically driven first 
oscillator [T] and phase equation for the autonomous second oscillator. If the 
external force is absent one can introduce relative phase of the oscillators 
9 = ipi — ifj2- In this case we obtain the classical Adler equation for two 

dissipative coupled oscillators 9 = — — — yusin(6') [IJ. This equation describes 

phase locking for the relative phase of the oscillators at | A| < 2/i and beats at 
|A| > 2/i. Phase locking of the oscillators in the case when external driving 
is absent occurs at there central frequency. 



3 Complete frequency locking of oscillators 
by the external force 

Now let us discuss the parameter plane (fi, 6) for the system First we 
shall find the domain of complete synchronization that corresponds to the 
exact phase locking of the both oscillators {ipi = "02 = 0). Substituting the 
last condition into Eqs.([5]) and excluding sin('?/'i — -02) in the first equation 
we obtain: 

2Q = 6sin(-?/'i), 
2fi-^ = /isin(V'i-^2). 
The system defines the boundaries of the synchronization region: 

b = ±2n, (7) 

Condition ([7]) defines the classical synchronization tongue in the plane(f2, b) 
(see FigUla). The cusp is at the point f2 = which corresponds to the fre- 
quency of phase locking of autonomous oscillators. However not all of the 
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regimes inside this synchronization domain are stable. Condition (|H]) sup- 
plements condition ([7]) and determines the width of a region of the stable 

synchronous regime. The center of this region coincides with the intrinsic 

A 

frequency of the second oscillator Q2 = And it's width is equal to a 
coupling constant /j, (Fig{T]a). 




Figure 1: Regions of phase locking (red color) and quasi-periodic dynamics 
(white color) for the driven system in the case of: a) phase locking of 
the autonomous coupled oscillators; b) quasi-periodic dynamics of the au- 
tonomous coupled oscillators. Q2 is the intrinsic frequency of the second os- 
cillator; be is the amplitude threshold of synchronization of the quasi-periodic 
oscillations; /i is the coupling parameter. 

The boundaries of the region of complete synchronization by the external 
force are given by the combination of the relations and ([H]). As a result, 
we obtain two types of the synchronization regions presented in Fig{TJa and 
FigJTJb. From relation ([7]) and (|8]) it follows that the first type of synchro- 
nization region is realized at |A| < 2/i, and the second type at |A| > 2/i. 
Hence, the first type of synchronization region corresponds to the case when 
the oscillators in the absence of the external forcing are phase locked, whereas 
the second type corresponds to the case of quasi-periodic dynamics. 

Comparison of FigJTJa and Fig{Tlb reveals some typical features of the 
synchronization of oscillator in the regime of beats. In this case an amplitude 
threshold for the complete synchronization of the quasi-periodic oscillations 
arises. From relations and ([S]) we find it as: 

6e = y - 2^1. (9) 

In accordance with relation ([8]), the frequency of the external force should 

A 

be close to the intrinsic frequency of the second oscillator = — to obtain 
the complete synchronization of the quasi-periodic regime. 
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4 Phase locking regime in autonomous oscil- 
lators 



System ([5]) can demonstrate various types of dynamics differing by tlie pliase 
portrait on tlie (■0i,'?/'2) plane. Typical phase portraits are shown in FigJ21 
There may be stable equihbrium points (Fig|2Ja); various types of attracting 
invariant curves (Figl2lb,c,e); or trajectories covering densely the whole phase 
plane (FiglJld). These types of dynamical behavior correspond respectively 
to the regimes of regular oscillations, two- and three-frequency tori in the 
initial system. 




Figure 2: Phase portraits of the system ([5]) at /i = 0.3 and A = 0.2: a) 
the complete synchronization of both oscillators, Q = 0.1, b = 0.3; b) the 
partial synchronization of phase of the first oscillator, Q = 0.75, b = 0.9; 
c) the partial synchronization of the relative phase of oscillators Q = 0.75, 
b = 0.4; d) the three-frequency torus, fl = 0.75, b = 0.75; e) the resonant two- 
frequency torus with winding number w = -, Q = 0.75, b = 0.79. Attractors 
are indicated by the red color. 

We will use next method (method of the charts of Lyapunov's expo- 
nents) for determining the type of the regime of system We calculate 
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all Lyapunov's exponents of the system Ai, A2 at each grid point in 
the parameter plane [Q, b). Then we color the points of the plane in accor- 
dance with values of the Lyapunov's exponents to visualize the domains of 
the corresponding regimes: 

1. P is the region of the stable equilibrium point (complete phase locking), 
Ai < 0, A2 < 0; 

2. T2 is the region of two-frequency quasi-periodic regime, Ai = 0, A2 < 0; 

3. T3 is the region of three-frequency quasi-periodic regime, Ai = 0, A2 = 
. 

Note that regime T2 (Figs|2jb,c,e) corresponds to an attractor in the form 

of two-frequency torus in the phase space of the original system. Accordingly, 

three-frequency torus is a phase portrait for regime T3 (Fig|2Jd). 

P 

The following number w = - is useful to classify the quasiperiodic regimes. 

Q 

Here p and q are numbers of crossings of the phase trajectory in steady regime 
with the sides of phase square {0 < ipi < 27r, < ip2 < Stt). Number w will 
be rational number (for which p and q are integer) for resonant two-frequency 
tori. In turn, for the three-frequency tori w is irrational 0. 

In Figjnia we present the charts of Lyapunov's exponents for the system 
([5]). It is plotted in the case when the autonomous oscillators are phase 
locked. Domain of the periodical regimes P presented in FigJSla corresponds 
to the analytical investigation of the complete synchronization presented 
above. System (EJ has three unstable and one stable equilibrium points 
in this region. The stable point defines the regime of the phase locking of 
coupled oscillators (Figj2ja). 

Two codimension-2 bifurcation point are disposed in (f2, b) plane (FigJSla). 
Their coordinates are 

The domains of two-frequency tori, three-frequency tori and periodical 
regimes adjoin at these points E|. These points determine the threshold of 



^Checking the zero Lyapunov exponents was carried out with an accuracy of 10^'^. The 
accuracy of calculating the main Lyapunov exponents was at least 10"''. 

^The value of w can be interpreted as a kind of winding number. Indeed, closing the 
region (0 < ■01 < Stt, < 7/12 < 27r) for both phase variables, we can obtain a "torus of 
phases" , for which w is a winding number. 

■^From the point of theory of bifurcations two lines of saddle-node bifurcations for 
stable points (O , ^ and lines of saddle- node bifurcations for stable and unstable invariant 
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I T2- two-frequency toms. 
T5- tluee-frequency toms. 

Figure 3: Charts of Lyapunov's exponents for the phase system ([5]) for next 
cases! a ) the phase locking of the autonomous oscillators, = 0.3 and A = 
0.2; b) the quasi-periodical behavior of the autonomous oscillators, /i = 0.3 
and A = 1.6. 

the amplitude of the signal at which three-frequency regimes are possible. 
Regions of three-frequency tori emerging from these points divide the 
domain of two-frequency tori into two typical regions: 

1. Partial synchronization of the first oscillator by external force. Winding 
number is = In this case the phase of the first oscillator -01 is 

fluctuating around some fixed value. While the phase of the second 
oscillator -02 demonstrates unlimited grow (Figj2jb). This regime is 
observed at large values of the parameter b. 

2. Partial synchronization of the relative phase of oscillators. Winding 

number is w = -. In this case the phases of both oscillators demon- 

1 ^ 
strate unlimited grow. However, relative phase of oscillators 6 = 1^1—1(^2 

is fluctuating around some fixed value (Figj2jc). This regime is observed 

at small values of the parameter b. 

curves (tori) are merging at the mentioned points [10,16-18,26]. Similar points for coupled 
rotation maps in Ref.[10] were called "saddle- node fan", due to the typical arrangement 
of resonant two-frequency tori of different order, see fig. 4. 21 in Ref.[10]. 
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The nature of these regimes is as follows. Small external force cannot 
destroy the strong phase locking of autonomous oscillators but it can syn- 
chronize both oscillators as a whole. Large force can destroy synchronization 
between the subsystems and then the first oscillator may be locked by the 
external force. Note, that partial synchronization of second oscillator is not 
possible. Second oscillator can only be synchronized by external force simul- 
taneously with the first oscillator within the domain P. 

In Figjnia one can see very narrow regions inside the T3 domain. Two- 
frequency tori are observed in these narrow regions. These are the high order 
resonant 2D tori situated on the surface of the corresponding three-frequency 
torus. The phase portrait (Figj2le) corresponds to the widest synchronization 

region. One can see that it is two-frequency torus with winding number 
1 



5 Quasi-periodic regime in autonomous oscil- 
lators 

Let us increase detuning parameter A and pass to the quasi-periodical regime 
in autonomous oscillators. Corresponding chart of lyapunov's exponents is 
shown in FigJSlb. The domain P of complete synchronization of both oscil- 
lators by the external force corresponds to that shown in Fig{Tlb. Complete 
synchronization regime demonstrates the threshold behavior with respect to 
the external force amplitude. The value of this threshold is determined by 
the formula (Q. 

The region of three-frequency tori transforms essentially. The amplitude 
threshold for 3D tori disappeared and they occupy the lower part of the 
parameter plane {Q,b). 

Region of the partial phase locking of the first oscillator w = - now looks 

like the traditional synchronization tongue locating within the domain of the 

three-frequency tori. The cusp of this synchronization region is at the point 
A 

= — — which corresponds to the intrinsic frequency of the first oscillator. 

Distinctive feature of this case is the possibility of existence of two- 
frequency torus with winding number w = -. It corresponds to partial 
synchronization of second oscillator. Domain of this regime has a form of the 

synchronization tongue whose top adjoins the region of the exact phase lock- 

A 

ing P while the cusp is at point = — which corresponding to the intrinsic 
frequency of the second oscillator. Thus, in the case under consideration, the 
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partial synchronization of the second oscillator which is not directly excited 
by the external force is possible. It constitutes important difference from the 
case when the autonomous oscillators are phase locked. In the latter case the 
coupling between oscillators is strong. And if external force has frequency 
Q ^ Q2 and its amplitude is small then only the relative phase can be locked. 
But if the frequency detuning A is large (the autonomous system of coupled 
oscillators demonstrated the quasi-periodical regime) then the second oscil- 
lator is not so strongly influenced by the first one. In that case external force 
can lock the second oscillator, while the phase of first oscillator is drifting. 

Regime of partial locking of the relative phase of the oscillators w = ^ 

at f2 > is mostly separated from the frequency axes by the wide region of 
three-frequency tori. 



6 Chart of two-frequency tori 

In FigJ3] one can see a set of narrow synchronization regions corresponding 

to different types of two-frequency tori. We use next method for its detail 

investigation. Let us fix a point in the parameter plane b). Then we solve 

numerically the system (|5]) for these parameter values, and find the number 

of intersections of p and q of the phase trajectory with the bottom and left 

sides of the phase square (0 < -01 < 27r,0 < 1^2 < 27r). (It is necessary to 

exclude irrelevant crossings associated with the possible oscillatory nature of 

the invariant curve.) Then we color the point in the parameter plane in a 

certain color depending on the type of the two-frequency torus which we can 

P 

characterize by the value w = -. Three- frequency tori are determined as 

Q 

non-periodic regimes. To obtain a chart we perform a similar procedure at 
each point of the grid in the parameters plane. 

In FiglHwe present chart of tori plotted using a method described above. 
This chart corresponds to the beats in the autonomous oscillators. In this 

case the widest synchronization region [w = -) corresponds to the partial 

A 

phase locking of the first oscillator. Its cusp is at the point at i7i = — — 

which corresponding to the intrinsic frequency of the first oscillator. The 

domain of partial phase locking of the second oscillator {w = -) looks like 

a continuation of the region of complete phase locking P. Its cusp is at 
A 

the point = — which corresponds to intrinsic frequency of the second 

oscillator. In turn, the point = ^2 is the center of whole system of the 
regions of two-frequency tori of higher order. 



11 




1110 1 1 



Figure 4: Charts of the tori for the case when the autonomous system demon- 
strates regimes of beats, fi = 0.3 and A = 1.6. Numbers indicate the values 
of winding number for the main two-frequency tori. 

7 Dynamics of the original system 

Phase equations ([5j) considered above are the approximation for the original 
system ([1]). They were obtained for the case of A << 1. Hence, it is important 
to clarify organization of the "frequency - amplitude" parameter plane for 
system ([T]) in the case when the value of the control parameter A is not small. 
Charts of Lyapunov's exponents for system ([T]) are shown in FiglSl They are 
plotted for A = 1. FiglSla corresponds to the case when the autonomous 
oscillators are phase locked. And FigjSlb corresponds to the case when the 
autonomous system demonstrates beats. 

Despite the large value of the control parameter A = 1, the picture has 
much in common with the case of the phase equations approximations. First, 
it concerns the configuration of the region of complete synchronization P. It 
is only slightly distorted as compared with FigJSl In the case when au- 
tonomous oscillators are phase locked one can see co dimension-two point in 
the region < (FigJSla). Clearly visible is a characteristic pattern of 
two- and three-frequency tori in its vicinity. The corresponding point in the 
high-frequency domain {Q > 0) which can be seen in Figl31a disappears to- 
gether with the adjacent region of three-frequency tori. Two-frequency tori 
and the system of very narrow regions of periodical regimes of higher order 
are observed in its place. These periodical regimes are absent in the phase 
equations. Periods of some regimes are indicated in FigJSlby numbers. These 
periods are defined numerically as the periods of cycles in the stroboscopic 



12 



-0.5 



Q 



0.5 



-0.5 



0.5 



a) 



Figure 5: Charts of Lyapunov's exponents for the periodically driven system 
of coupled van der Pol oscillators ([I]): a) phase locking of autonomous oscil- 
lators, = 0.1 and A = 0.1; b) beats of autonomous oscillators, /i = 0.1 and 
A = 0.6. 

section of the system ([T]). 

A chart plotted for the case of beats in autonomous oscillators (Fig|5lb) 
also largely corresponds to the case of the phase equations (FigJSlb). First, we 
should note the presence of amplitude threshold for the complete synchro- 
nization of quasi-periodic oscillations. Also in the lower part of the chart 
one can see the region of three-frequency tori, inside of which the regions of 
different resonant two-frequency tori described in the previous section are ob- 
served. However, the high-frequency codimension-two point also disappears. 

Another new effect in comparison with the case of phase equations is an 
appearance at FigjSlb of a set of additional "fans" of the regions of two- 
frequency tori on the right side of the chart. They are not associated with 
the region of complete synchronization. These types of two-frequency tori 
are not observed in the phase equations. 

Stroboscopic map is a convenient way to visualize the type of the dynam- 
ical regime. It is a shift map along the trajectories of the periodically driven 
system taken over the period of external force. The result is an orbit - a set 
of points embedded in the space of dynamical variables of the system ([T]) (x, 
X, y, y). For their visualization we use projections on the plane of the dy- 
namical variables of the first and second oscillators (x, x) and (?/, y). These 
illustrations are given in the first two columns in FigjG] The third column 
shows the same projection on the velocities plane of the two oscillators (i, 
ij.) It is a kind of Lissajou's figures. For the orbits from the first two columns 
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of FigE] let us determine the phases of oscillators tpi and ip2, as the angles 
at which the points of the orbit are visible from the origin. The trajecto- 
ries on the phase plane {tpi, v?2) obtained using this method are presented 
in the fourth column in FiglHl FigJHla illustrates the three-frequency torus 
in the system ([1]). The trajectory in this case covers densely the plane of 
the oscillators phases {tpi, ip2)- FigElb illustrates a resonant two-frequency 

torus with winding number w = -. This two- frequency torus arises from a 

three-frequency torus with a small change in the amplitude of the external 
force. 




Figure 6: Stroboscopic section of the forced system of coupled van der Pol 
oscillators ([T]) in the projection on the plane of the variables of the first 
and second oscillators, in the form of Lissajous figures, and portraits of the 
plane of the oscillators phases in the stroboscopic section. Parameters are a) 
B = 0.45, h) B = 0.5; A = 1, A = 0.1, n = 0.1, n = -0.25. 



8 Conclusion 

It is necessary to distinguish the cases of phase locking and quasi-periodic 
dynamics with incommensurate frequencies in autonomous oscillators when 
we investigate synchronization of two-frequency oscillations by the external 
force. In the first case regimes of complete synchronization of two oscillators, 
partial synchronization of relative phase of oscillators and partial synchro- 
nization of first oscillator are prevailing. Appearing of three-frequency tori is 
possible only above some external forcing amplitude threshold. Simultaneous 
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exact locking of phases of both oscillators by harmonic external force is also 
possible in the case of quasi-periodic dynamics of autonomous oscillators. 
However, now this effect demonstrates amplitude threshold. There exists 
domain of regime of partial synchronization of second oscillator as amplitude 
decrease: it has a form of synchronization tongue in the parameter plane of 
frequency versus amplitude which is surrounded by system of synchroniza- 
tion regions of resonant two-frequency tori of high order. The amplitude 
threshold for three-frequency tori disappears in that case. The method of 
charts of dynamical regimes is an effective method for the analysis of such 
systems. With this method we color the parameter plane in different colors 
in accordance with winding number, which can be attributed to each reso- 
nant two-frequency torus. This method may be used for the further analysis 
of other similar systems. For increasing values of the control parameter of 
the original system of coupled van der Pol equations the main features of the 
organization of the parameters plane are preserved. In particular these are 
the features characteristic for the transition from mode locking to the beats 
in autonomous oscillators. However, at large amplitudes the three-frequency 
tori disappear. In their place the domain of two-frequency tori along with 
system of periodic regimes comes. The latter regimes have no analogue in 
the phase equations. 

We thank S.P. Kuznetsov for help and discussions. This research was 
supported by the grant RFBR No. 09-02-00426. 
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